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We calculate the s-wave part of the pion-nucleus optical potential using a unitarized chiral ap- 
proach that has been previously used to simultaneously describe pionic hydrogen and deuterium 
data as well as low energy irN scattering in the vacuum. This energy dependent model allows for 
additional isoscalar parts in the potential from multiple rescattering. We consider Pauli blocking and 
pion polarization in an asymmetric nuclear matter environment. Also, higher order corrections of 
the irN amplitude are included. The model can accommodate the repulsion required by phenomeno- 
logical fits, though the theoretical uncertainties are bigger than previously thought. At the same 
time, we also find an enhancement of the isovector part compatible with empirical determinations. 
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I. INTRODUCTION 

The problem of the missing repulsion in pionic atoms has attracted much attention in the past [1-10] and recently 
[11-17] and was further motivated by the discovery of deeply bound pionic atoms at GSI [18-21]. 

Due to the repulsion of the s-state pion in nuclear matter, the tt" wave function is strongly repelled and overlaps 
only little with the nucleus. The wave function tests mainly the peripheral zone of the nucleus and, thus, nuclear 
matter at less than nuclear density. However, even at half the nuclear matter density difficulties in the theoretical 
description persist. From phenomenological fits to pionic atoms reaching from C to Pb, a strong repulsion is needed 
for a consistent description of the combined data. However, theoretical calculations consistently failed to deliver this 
"missing repulsion" (see, e.g., [5]) although there has been recent progress [14]. 

The s-wave pion-nucleus optical potential is the basic input for a calculation of the s-levels of pionic atoms. 
Usually, the optical potential is calculated for infinite nuclear matter as a function of the Fermi momentum. Explicit 
calculations for finite nuclei done in Ref. [22, 23] provide a prescription to pass from nuclear matter to finite nuclei: 
the s-wave part of the potential is provided by the corresponding nuclear matter results changing p to p{r) (local 
density approximation), while for the p-wave the prescription is slightly more complicated. 

The s-wave pion optical potential 2wl4pt(r) = Hs(r) is closely connected to the s-wave pion sclfencrgy which is 
usually [1] parametrized as 



Ils{r) = -47r 



boiPp + Pn) + ( 1 + ) bliPn - Pp) + ( 1 + ^ j Bo{p){pp + Pnf 



(1) 



where the density p is a function of the radial distance, p = p{r), given by the density profile of the nucleus. From 
this expression the sensitivity of the selfenergy to the isoscalar bo becomes visible, as in symmetric nuclear matter 
the isovector term bi vanishes. However, heavy nuclei such as |2*Pb recently used in experiments [19-21] contain 
more neutrons than protons; it is therefore interesting to study asymmetric matter, in particular with respect to a 
possible renormalization of the isovector bi [15, 24-26]. The last term in Eq. (1) takes into account corrections from 
higher order in density. This quantity has also an imaginary part due to pion absorption, which is mainly a two-body 
process, and the imaginary part of the optical potential determines the width of the pionic atom. 

Traditional fits to pionic atoms [22, 27, 28] provide the set of parameters displayed in Tab. I. Although the sets 
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TABLE I: Typical fits of pionic atom dat. 


3,. 


Ref. 


bo [m;^^] 


bi [m-i] 


Bo [m-4] 


[27] 


-0.0045 


-0.0873 


-0.049 + i 0.046 


[28] 


-0.0325 


-0.0947 


0.002 + i 0.047 


[22] 


-0.0183 


-0.105 


i 0.0434 



of parameters are quite different from each other they result in similar pion self energies at p = po/2, half the nuclear 
density. Therefore these sets are not contradictory but tell us that the pionic atom data require this magnitude of 
selfenergy at po/2. This equivalence of pion optical potentials using the concept of peff = /fo/2 was early established 
in [29, 30]. Furthermore, Tab. I suggests that the smaller value of |6o| in Ref. [27] needs to be compensated by a 
large negative real part of the p^-term Bq; thus, corrections of higher order in the density are important. 

The model of Ref. [31] is of interest in this context as a good part of the ttN vacuum isoscalar is generated by the 
multiple rescattering of the dominant Weinberg- Tomozawa term of the ttN interaction of isovector character. This 
realization is important because rescattering terms are appreciably modified in the nuclear medium. Indeed, the Pauli 
blocking in the intermediate nucleon states is well known to generate a repulsion, the Ericson-Ericson Pauli corrected 
rescattering term [1]. On the other hand, the pion polarization due to particle-hole (ph) and A- hole (Ah) excitation 
of the intermediate pions also produces corrections and accounts for the imaginary part of the potential from pion 
absorption [5, 32]. 

Another point is the energy dependence of the TriV interaction [14]. Ref. [31] focuses on the precise determination of 
the scattering lengths but also provides the energy dependence close to threshold. For pionic atoms where the pion is 
practically at rest with respect to the nucleus this is still relevant due to the Fermi motion of the nucleons. Note that 
in this context, the vacuum model [31] already contains certain information about the nucleon-nucleon correlations as 
one of the fitted data points has been the 7r~-deuteron scattering length. The deuteron wave function that enters the 
theoretical description provides the NN momentum distribution and allows for an inclusion of the Fermi motion in 
the deuteron. The issue of the energy dependence is a relevant one and in the medium it induces corrections which, 
due to the smallncss of the bo parameter, have an effect similar to a renormalization of bi [14]. 

On the other hand there are some medium corrections coming from vertex corrections, off-shell effects, and wave 
function renormalization which, if desired, can also be recast as renormalization of bi and bo. We shall also introduce 
novel terms in the pion selfenergy related to the TV* (1440) decay into Ntttt, with the two pions in a scalar isoscalar 
state. This mechanism has already been used in [33] to estimate some uncertainties in the study of the 7r-deuteron 
interaction. 

Another novelty in the present work is that we shall start from a free model for ttN scattering which is constructed 
using a chiral unitary approach, incorporating the lowest order (LO), and the needed next to lowest order (NLO) 
chiral Lagrangians, together with multiple scattering of the pions [31]. 

The vacuum model from Ref. [31] will be modified in various steps in the medium: In Sec. II A Pauli blocking of 
the intermediate nucleonic states, together with the appropriate spectral function for the intermediate pions, will lead 
to non-linear corrections in the density with preliminary numerical results given in Sec. III. Also, a self consistent 
calculation is presented in Sec. Ill A where the overall pion s-wave selfenergy serves as an input for the intermediate 
pions in the ttN loops. In Sees. IV, V, the diagrammatic model will be extended to the above mentioned higher order 
vertex corrections. Final numerical results are provided in Sec. VII. 

II. LOW ENERGY PION NUCLEON INTERACTION IN VACUUM AND MATTER 

The vacuum ttN isoscalar term bo is around ten times smaller than the vacuum isovector &i-term and its precise 
determination is a complex task due to large cancellations in the amplitude. With the advent of new experimental data 
[34-38] for the iT~p tt~p, Tr~p Tr'^p, and ir^d ir^d scattering lengths from pionic hydrogen and deuterium, 
theoretical efforts in several directions have been made to precisely determine the parameters of low energy TriV 
scattering. In this context, 7r~-dcuteron scattering at threshold plays an important role as the complex scattering 
length a.^-d puts tight constraints on the size of bo- 

Pion deuteron scattering has been recently treated in chiral perturbative approaches [33, 39] including also correc- 
tions from isospin breaking [40] and effects like Fermi motion [41]. These and other higher order corrections have been 
taken into account in another theoretical framework in Ref. [42]. In the extraction of the strong scattering lengths 
from experiment, special attention has to be paid to the Coulomb corrections in the extraction of the scattering 
lengths from pionic hydrogen [43-45]. 

In the present study we rely upon a recent study on low energy ttN scattering in s-wave [31] which is summarized 
below. This model simultaneously describes the available data at threshold from pionic hydrogen and deuterium and 
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FIG. 1: Rescattering of the N system generated by the Bethe-Salpeter equation. 

also low energy ttN scattering. In a restriction to the coupled channels Tr^p, n'^n, and n^n the ttN s-wave amplitude 
T(^/s) is unitarized by the use of the Bethe-Salpeter equation 

r(Vi) = [1 - V{V~s)G{V~s))] ViyTs). (2) 

Here, the kernel V is given by the elementary isovector interaction from the Weinberg- Tomozawa term of the LO 
chiral Lagrangian [46-48], 



iV-s) = -a, jj, {2V-S - - M,) ^ ^ . (3) 

The ttN loop function G in Eq. (2) provides the unitarity cut and is regularized in dimensional regularization with 
one free parameter, the subtraction constant Ut^n [31]. In Eq. (3) the coefficients provide the transition strength 
of the coupled channels i,j [31] and Mij, Eij are the nucleon masses and energies. In Fig. 1 we show a diagrammatic 
representation of the BSE equation (2), including also the ttttTV channel which is also incorporated in [31]. 

In the framework of the heavy baryon approach the vertices are factorized on-shell, see Eq. (3), because the off-shell 
part of the vertices in the loops can be absorbed renormalizing the lowest order tree level amplitude [49]. However, 
we will see in Sec. IV that in a nuclear matter environment these renormalizations are modified leading to finite, 
density dependent corrections of the amplitude. 

The multiple rescattering which is provided by Eq. (2) generates isoscalar pieces from the isovector interaction 
providing a large bg term. However, it is known [50-52] that the NLO chiral Lagrangian is a necessary ingredient in 
ttN scattering at low energies. In order to provide the necessary degrees of freedom in the model, the isoscalar s-wave 
piece with the chiral coefficients Ci in the notation of Ref. [52] , 

v., - v., + % ^ - 2c, j 2,,^ , (4) 

is added to the kernel of Eq. (2). The term C3 in ref. [52] has been taken as C3 to^, consistently with the approach 
of Refs. [53, 54] which uses the on-shell values for the vertices in the scattering equations. The free fit parameters 
up to this point are the subtraction constant aj^N and the two combinations of q from Eq. (4) , as well as a damping 
factor parametrized with f3 as discussed in Ref. [31] which is of no relevance here because we stay close to threshold. 
There are further refinements of the model, described in detail in Ref. [31], such as the inclusion of the ttttN two-loop 
diagram which introduces one additional fit parameter, 7, from the small real part of this loop. 

In order to include the complex pion-deuteron scattering length a^r-d in the data fit, one has to employ the 
elementary ttTV scattering model described above in the framework of a three body process. In Ref. [31] this has been 
carried out by using the ttN amplitudes in a Faddeev multiple scattering approach. The interesting point is that the 
impulse approximation vanishes making the double rescattering off the two nucleons the dominant term. This term 
is sensitive to the isoscalar amplitude so that the experimental scattering length a^^-d provides valuable information 
on the vacuum 60 term and sets tight constraints on it. 

Additional corrections of higher order Ref. [42] in nd scattering such as absorption, dispersion, the influence of the 
A(1232), and Fermi motion have been treated in a separate Feynman diagrammatic approach, together with other 
corrections from the literature, see Ref. [42] and references therein. Once these various corrections are included in 
and, the model parameters are fixed from data, namely the scattering lengths aT^-p^T^-p, o,Ti-p^Ti'>m o-Ti-di a-nd low 
energy ttN data from [56]. The parameter values are quoted in the left column of Tab. II. The values of Ci from Eq. 
(4) are in agreement with other works [52]; furthermore, the isospin violations found in the study qualitatively agree 
with Ref. [57]. In the following, we concentrate on the in-medium modifications of the approach. 

A. The model in nuclear matter 

The s-wave ttN — > ttN vacuum model from Ref. [31], summarized in Sec. II, provides the driving interaction of the 
7r~ with the nucleus. In order to obtain the pion selfenergy H^ from Eq. (1) of the tt^ in asymmetric nuclear matter 



FIG. 2: Diagrammatic representation of the vr selfenergy from s-wave interaction with nucleons. 

with proton and neutron densities Pp and Pn (k^p, kp the respective Fermi momenta), the tt^N —>■ N ampUtudc T 
is summed over the nucleons in the Fermi sea as schematicahy indicated in Fig. 2. The s-wave selfenergy for a tt^ at 
momentum (fc°,k) with respect to the nuclear matter rest frame reads 



(5) 



where Pp^„ are the nucleon momenta. Due to the breaking of Lorcntz invariancc, the amplitudes Tj^-p,n depend 
independently on the components of (P°,P), the total 4-momentum of the TriV system in the nuclear matter frame, 
namely P° = + Ep^n{Pp,n) and P = k + pp_„. The factors of 2 in Eq. (5) account for the sum over the nucleon spins. 
Note that Eq. (5) allows for isospin breaking by using different masses for particles of the same isospin multiplet. 
In analogy to the vacuum case, T^-p and r^-„ are given by the solutions of Bethe-Salpeter equations (BSE) 

T(P", P; p)=[l- V{Vs)G{P\ P; p))] V{V^) (6) 

where s = (P°)^ — P^ and the loop function G is modified as described below. In Sec. IV we will apply in-medium 
changes also to the kernel V from off-shell parts of the vertices and other sources. For the charge C = sector, the 
BSE is represented by (2 x 2) matrices accounting for the coupled channels 7r~p and Tr'^n. For the 7r~n interaction 
there is only one channel. 

The diagonal matrix G from Eq. (6) contains the loop functions Gt^n which have been formulated in dimensional 
regularization in Rcf. [31, 54] for the vacuum case. Alternatively, one can use a cut-off scheme [54] with A the three 
momentum cut-off. The vacuum G^at is then given by 

r (pO-p\_ r d'^q Mn 1 1 

with a cut-off for the three- momentum integration A = 1 GeV and m (Afjv) being the 7r~,7r° (p, ?^) masses. Over 
wide energy ranges, a change in A can be written as an additive constant to the real part of G^^n- Therefore, we 
have denoted a separate piece Otttv in Eq- (7) in the same way as in Rcf. [54]. For the free case the propagator in 
the cut-off scheme agrees with the propagator from dimensional regularization over a wide energy range by choosing 
the appropriate subtraction constant. In the nuclear medium with Lorcntz covariancc explicitly broken, a cut-off 
scheme is more convenient in order to implement the in- medium dressing. Thus, we will employ the propagator from 
Eq. (7) in this work. This requires a refit of the vacuum data. The values of the model parameters with the cut-off 
propagator from Eq. (7) instead of dimensional regularization are displayed in Tab. II on the right hand side. The 
new fit shows that the model is insensitive to the used regularization scheme. Parameters, x^, and predictions for 
isoscalar and isovector terms &o and bi are stable. For notation of the parameters, see Sec. II. In Tab. II, ajrAf is the 
subtraction constant of the loop in dimensional regularization and Atttv the subtraction constant from Eq. (7). The 
more important parameters are the Ci and a7r7v(a7r7v)- The real part of the ttttN loop (7) is tiny at threshold. For 
pionic atoms, the damping factor /3 of [31] which is more important for the higher energy nN data is of no relevance 
because the cm. energy of nN due to Fermi motion in the nucleus is small. 

The two major medium modifications of Gj^n are the Pauli blocking of the nucleon propagator and the polarization 
of the pion. The corresponding diagram is displayed in Fig. 3. For the amplitude of the in-medium nN loop function 
a similar expression as in Ref. [58] is obtained. Here, we give the generalization to asymmetric nuclear matter for the 
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TABLE II: Global fits to pionic hydrogen, deuteron, and low energy vrA'^ scattering data, using dimensional regularization from 
Ref. [31] and cut-off scheme. Also, the resulting feo,&i are shown. 





Til TTl f^CF 
I #111111 




fitted data (Vi) 


1104-1253 MeV + threshold 


1104-1253 MeV + threshold 


2 
Xr 


51/(2 - 10 + 3) 2.2 


48/(2 • 10 + 3) ~ 2.1 


r 1 

QttJV [-J 


1 1 yi Q _i_ n 1 no 




O-ttN [MeVJ 




o nott _i_ 1 OQ 


2ci - C3 [GeV"^] 


— 1.539 ± 0.20 


— 1.487 =b 0.20 


C2 [CjeV \ 


— 2. DO ( ± y).2Z 


— 2.000 ± U.2z 


p [MeV J 


n nnoT/i i _i_ i c i n— 4 
U.UUz/4i zb i.O - iU 


A AAOTttO _l_ 1 C 1 A~4 


7 [10-^ • ml] 


5.53 ±7.7 


6.27-7.8 


\ (.^7V~ p — >7r~f)) 


3 


3 




< 1 


< 1 


X^ia^-d) 


8 


7 


bo [10-* m-':] 


-28 ± 40 


-29 


&i [10-* m;^] 


-881 ± 48 


-883 



\ . hole N 

X 1^^-^ 

FIG. 3: In-medium correction of s-wave ttA*' scattering: Renormalization of the pion and Pauli blocking of the nucleon, symbol- 
ized by a crossed propagator. The pion p-wave selfenergy stands for resummed p/i, A/i insertions and includes A^'A*', AA, AA 
short-range correlations. 



TT p, Tr'^n, and tt n loops. With N ^ p,n and tt^ = tt , tt" 



,max 

0(fc#-|P-q|) 



|qcm|)- 



g(|P-q|-fcf) 

pn - - £;Ar(P - q) + ie ^ P" - 9" - Pjv(P - q) - «e 



M 



AT 



du) 



2uj 



- 5^.(tj,q;pp,p„). (8) 



The cut-off in the vacuum model is applied in the t:N cm. frame, as required by the vacuum model, whereas Eq. (8) 
is defined in the nuclear matter rest frame. Since in the free case gmax is given in the cm. frame we boost q to this 
frame and demand it to be smaller in modulus than q^^- We have 



P" 



1 



P q 



q" 



(9) 



where s = (P°)^ — P^. In Eq. (8) we have also taken into account the hole part of the nucleon propagator as in Ref. 
[55] which can play a role at the low pion energies we are studying. This term has been neglected in Ref. [58] which is 
justified at higher energies. The pion spectral function 5'.^.. is different for Tr~ and tt^ for asymmetric nuclear matter. 
For S we include the particle hole (ph) excitation and NN short-range correlations as described in the next section. 

In the model from Ref. [31], the A(1232) has been explicitly taken into account in pion-deuteron scattering, leading 
to corrections in the ird scattering length which by itself sets constraints on the vacuum isoscalar amplitude. In the 
present situation we can take the corresponding effect into account by including also the A-hole (Ah) excitation in 
the pion selfenergy; in fact, closing the nucleon lines of the deuteron in the A-box and A-crossed box diagrams of 
Ref. [31] one obtains a pion selfenergy corresponding to Fig. 3 substituting the ph by a Ah excitation of the pion. 
The A''* (1440) Roper-hole excitation can be in principle also included in the pion selfenergy but has been found small 
in Ref. [58] for low energy pions. However, in Sec. VII B the Roper resonance will be included in a different context 
based on the coupling of the Roper to a scalar-isoscalar pion pair. 



FIG. 4: Integration over the Fermi sea of the medium -kN amplitude. The crosses represent Pauli blocking of the nucleon 
propagators and the large dots, the p-wave pion selfenergy. 



In Ref. [58] Pauli blocking for the intermediate ttttN loop - see Fig. 1 - has been included for the imaginary 
part. In the present case the pion has very little momentum in the 'kN cm. frame and the system is below the ttttN 
threshold, hence, even in vacuum the imaginary part of this term is zero and thus no change is required. Hence, the 
imaginary part of the amplitude from the ttttTV intermediate state is zero in our case and the contribution to the real 
part of the amplitude is in any case negligible. 

Combining all the ingredients of the in-medium model, the s-wave pion selfenergy can be symbolized by the diagram 
in Fig. 4: The in- medium propagator from Eq. (8) is used in the BSE (6), and the remaining integral over the Fermi 
seas from Eq. (5) corresponds to closing the nucleon line. 



B. Pion polarization in asymmetric nuclear matter 

The spectral function of the pion tt^ (tt^, tt^, tt^) at momentum [q^ , q) from Eq. (8) is given by the imaginary part 
of the propagator, 

For the pion selfenergy inside loops the p-wavc part is dominant because q is a running variable and 117^; a q^ . The 
s-wave part will be included in the self consistent treatment in Sec. Ill A. For the selfenergy we take into consideration 
the {ph) — (ph) short-range repulsion parametrized in terms of the Migdal parameter, which is chosen g' = 0.7, 

TT ^ (^ + ^\^ j^2t^ 2 t^7r,(g°,q;Pp,Pn) 

n7r,(g ,q;pp,p«) = \^-, — ^ (g)q ; ■ (n) 



For a diagrammatic representation of the pion selfenergy, see, e.g., Ref. [9]. The Lindhard functions for asymmetric 
matter for [ph) and (Aft.) excitations, evaluated below, are added in Eq. (11), U = [/(p'') -|- J/^^''). Note that we have 
here for simplicity assigned the same g' to {ph) and (Aft) excitations. For the form factor that takes into account the 
off-shell pions coupling to ph or Aft we have chosen the same function F{q) = A^/ (A^ + q^) with A = 0.9 GeV. 

The Lindhard function for symmetric nuclear matter, U{q,kp), can be found in the literature, e.g. in Ref. [59], 
and here, we concentrate on an extension to asymmetric matter (see also [60]). In the non-relativistic reduction, the 
Lindhard function for pions turns out to be 



C/..(9,4:4) = 4 J 



(27r)3 



e(fci,-iki)e(ik + qi-fc|,) , e(fc2 -iki)e(ik-qi-fci,) 



+ e(k) - e(k + q) + ii] -g" + e(k) - e(k - q) + ii] 



(12) 



The first term is the contribution of the forward going ph excitation (direct term) and the second term the pion 
crossed-term selfenergy. The index 1 (2) labels the Fermi sea corresponding to the hole (particle) part of the direct 
and the particle (hole) part of the crossed contribution. E.g., for a tt~ , kp = kp and kp = kp. For a 7r+, kp = kp, 
kp = kp. The integral (12) can be solved analytically. For this, we split the ordinary Lindhard function from Ref. 
[59] in direct and crossed part by ?7(g°, q, A:_f) = Ud{q^,(l, kp) +Uc{q^,q,kF) with Uc{q'^ , kp) = Ud{—q°,q^,kp) and 



where p = 2/(37r^)fc|, and A/jv the proton or neutron mass. Evaluating the integral in Eq. (12) one obtains for the 
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FIG. 5: Real and imaginary part of the pion propagator for symmetric and asymmetric nuclear matter at a pion momentum 
of 500 MeV. The position of the quasielastic pion peak in vacuum, at — 518 MeV, is indicated with the arrows. The 
asymmetric matter corresponds to the ratio of n to p in I^^Pb. 



ph Lindhard function in asymmetric matter 

C/d(g°,q,fcI^) + C/c(g°,q,fc^), 
C/d(9°,q, fc?^) + C/c('7°,q, k^p), 

^{U{q°,q,kP,) + U{q°,q,k'^p)). (14) 

These are the expressions to be used in Eq. (11). The result in Eqs. (12) and (14) is in agreement with Ref. [61], 
correcting a typographical error in their Eq. (A. 5). 

For the Ah Lindhard function U'-'^''\q, k^p, fc^,), no new calculation is required, as the A always plays the role of a 
particle and is not affected by the Fermi sea. It is therefore sufhcient to split U^{kp) from Ref. [59] into its charge 
states and direct plus crossed parts, and use as argument the kp that corresponds to the hole part, 

C/f''^g°,q;Pp,P«) = it/r)(g",q;fc^) + \ui^^\q\c^,ep) + \u^^^Hq',ci;k"p) + ^C/f'')(g°, q; fc^^), 

C/^^''('Z°,q;Pp,P«) - i(c/^^(g^q;fc?.) + L/i^'')(g^q;fc?,) + C/i^'')(g^q;fc^^) + C/f'')(gO^ (15) 

Analytic expressions for the direct and crossed part of the Ah Lindhard function can be found in Ref. [59] . 

In order to see the effects of asymmetric nuclear matter we plot the pion propagator for normal nuclear density po = 
0.483 which corresponds to kp = 268 MeV for symmetric matter. For asymmetric matter we set kp = 1.154 kp 
which corresponds to the ratio of neutron rich nuclei such as |2*Pb. Then, po = Pp + Pn is obtained with kp = 247 
MeV and kp = 286 MeV. In the plots in Fig. 5 the propagator from Eq. (10) for pions at |q| = 500 MeV is shown. 
The 7r° in asymmetric matter is very similar to the case of symmetric matter. The tt~ shows some minor deviations. 



UlT'iq',ci;pp,p,,) = 
C/^-''(9°,q;Pp,Pn) = 
ul'!,''\q°,q;Pp,Pn) = 



III. NUMERICAL RESULTS 



In Fig. 6, the real part of the s-wave pion sclfenergy from the full model and from several approximations is plotted. 
The solid lines show the results for the model from Sec. II A in symmetric and asymmetric nuclear matter, with 
and without the p-wave renormalization from Eq. (11) of the pion propagator in the intermediate ttA^ loops. For 
the cases with asymmetric matter, the x-axis is given by fc^. The neutron Fermi momentum is then chosen to be 
kp = 1.154fc^. This ratio corresponds to the ratio of neutron rich nuclei as |2*Pb with kp = 241 MeV and kp = 278 
MeV. The selfenergy in asymmetric nuclear matter is larger than in symmetric matter which can be easily understood 
from the large and positive term (— 47r)6i(p„ — pp) from Eq. (1). 

The effect of Pauli blocking in the intermediate loops of the s-wave rescattering (sec Eq. (8)) can be taken into 
account by the Ericson approximation [1, 5] 



(16) 
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FIG. 6: Real part of the s-wave pion selfenergy for the pion at rest. Note that for asymmetric nuclear matter, kp of the proton 
is plotted on the abscissa and we always take kp = 1.157 fc^. Fit results to pionic atom data from Refs. [22, 27, 28] are also 
plotted. Theoretical calculation of Ref. [5] indicated as "Garcia et al." (dashed line). 




FIG. 7: Imaginary part of the s-wave pion selfenergy for the pion at rest. Phenomenological fits as in Fig. 6. 
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As pointed out in Ref. [5] the quantities Ai^2 are related to the vacuum isoscalar and isovcctor bo, bi terms (generated 
from rescattering, not the elementary ones) for which we take from Ref. [31], 



&o, vac = -0.0028 m-'=- 



2Ai 



1 + 12^ m^' 



1 2A 

6i, vac = -0.0881 m-^= - ^ — i. (17) 



mjv 



With these values and Bq ~ (no pion medium modification) one obtains from Eq. (1) the dotted curve in Fig. 6 
upper left panel. Adding the approximate medium change of b^ from Eq. (16) according to feo = ^o, vac + ^bg, the 
dashed curve is obtained. Thus, the tp approximation is not sufficient, whereas the inclusion of Abg leads to a good 
agreement with the rescattering model. This shows also that effects from Pauli blocking in more than one loop in the 
ttN rescattering of the ttN amplitude are small, because Eq. (16) corresponds to exactly one Pauli blocked loop in 
the rescattering scries [5]. 

Next, we compare to asymmetric nuclear matter but still without pion modification. This is displayed in Fig. 6 
upper right panel. Now, the isovector term contributes and we can derive a similar approximation as Eq. (16) for the 
bi renormalization in nuclear matter, 

A6i(fc^) = (2A1A2 - Ai) (18) 

by simply comparing the isospin structure of TriV scattering at one loop. The result from Eq. (1) using 60 = &o. vac+A6o 
and bi ~ 5i, vac is indicated as the dotted line. As the dashed line we plot the result from Eq. (1) using bo = bo, vac+A6o 
and 61 = bi, vac + A61. Obviously, the correction from Eq. (18) is small. However, in Sec. IV we will find additional 
vertex corrections that will modify appreciably the isovector strength of TriV scattering. 

When including the pion renormalization in the model according to Eqs. (8,10) the real part of the s-wave pion 
selfenergy for symmetric and asymmetric nuclear matter decreases as shown in the two lower plots of Fig. 6. We can 
compare to Ref. [5]. For this, we take the final values for Bo from there. Bo = 0.032 + i 0.040 m^^. Note that this 
is only qualitative because we do not take the density dependence of Bo from Ref. [5] into account but use a mean 
value. The values from [5] for bo and bi are — 0.013 and — 0.092 to~^, respectively. With these values and adding 
A60 from Eq. (16) to bo, the selfenergy is calculated according to Eq. (1) and plotted in Fig. 6, lower left panel 
for symmetric nuclear matter. In the same plot the s-wave selfenergy from fits to the bulk of pionic atom data from 
Refs. [22, 27, 28] with the values given in Tab. I is shown. Both the present model and results from Ref. [5] are 
systematically below the phenomenological values. Neither the present model nor Ref. [5] reach the required size for 
the real part of Us and thus the problem of missing repulsion persists. 

The imaginary part of the pion s-wave selfenergy is displayed in Fig. 7. The result from Ref. [5] (dashed line) 
agrees well with the phenomenological values from Refs. [27] and [28] (gray band) whereas the present model shows 
a 30% discrepancy. 

The differences between the results from Ref. [5] and the present calculation (dashed vs. solid line for the symmetric 
matter case including the pion renormalization) should be attributed to a different input used in [5], such as form 
factors plus the fact that extra crossed terms of character (smaller than those incorporated here) were also evaluated 
in [5]. The larger repulsion from [5] can be partly explained by the large vacuum |6o|, |&i| used there, whereas nowadays 
values for bo, compatible with zero as in Eq. (17), are regarded as more realistic. 

A. Self consistent treatment of the amplitude 

For the pion polarization in intermediate ttN loops, so far only the p-wave pion selfenergy has been taken into 
account. For the s-wave part we can include the selfenergy determined in the last section in a self consistent approach. 
For this, the 7r~ selfenergy Us from Eq. (5) is included in the pion propagator from Eq. (10). Additionally, the 
selfenergy is resummed so that it can be included in the same way as the p-wave selfenergy IIp^ in the pion propagator, 

jj^ = - (19) 

(q0)2 -q2 -ml^ - Up,^{q° ,q; pp, pn) -Ils{q° = m^,q = 0;pp,p„)' 

We have approximated here the energy and momentum dependence of lis by the static case ((7° = ttitt, q = 0). Solving 
for lis by iteration one obtains the results in Tab. Ill for asymmetric matter. As in Sec. Ill we set kp = 1.157 kp 
and show the results for kp = 213 MeV and kp — 241 MeV which corresponds to densities of around po/2 and 
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TABLE III: Self consistent treatment of the s-wave selfenergy ns(g" — m^,q = 0) in [MeV^] for asymmetric matter. To the 
left, the case with kp = 213 MeV, to the right kp = 241 MeV. Three iteration steps are shown. 





Re(ns)[2l3McV] 


Im(ns)[2i3McV] 


Re(ns)[241 McV] 


Im(ns)[241 MoV] 


Step 


2470.4 


-570.8 


3423.6 


-1233.8 


Step 1 


2503.9 


-562.4 


3491.3 


-1207.4 


Step 2 


2504.3 


-562.1 


3492.2 


-1205.8 




FIG. 8: The vertex tadpole at 1//* (dia gram 1) and the corresponding medium diagram at (diagram 1'). 

Pq. Three iteration steps are shown with step being the selfenergy without iteration. Comparing the size of lis 
from Figs. 6 and 7 with from the propagator, the result is expected to change only little. Indeed, the iteration 
converges rapidly and changes are small. At this point one can improve the calculation by evaluating the s-wave pion 
self energy not in the approximation (g*^ = m-^, q = 0) as in Eq. (19), but with the full , q dependence: it is known, 
at least for the vacuum case, that the isoscalar ttN amplitude is small at threshold but then grows rapidly at finite 
scattering energies. Taking only the g^-dependence — the q dependence is small — the self consistent calculation 
delivers indeed a larger change than before, of about 10 % of additional repulsion at p = po/2. 

IV. HIGHER ORDER CORRECTIONS OF THE ISOVECTOR INTERACTION 

In this section additional corrections are introduced that go beyond the medium modifications from Sec. II A, 
namely medium corrections affecting the kernel of the Bethe-Salpeter equation itself. In our model the kernel is 
given by the Weinberg- Tomozawa isovector ttN ttN transition and the NLO isoscalar ttN — ^ ttA^ transition. 
Considering vertex corrections of the rescattering is advantageous because it allows to include higher order corrections 
to the Ericson-Ericson rescattering piece, that is a large source of isoscalar strength. The corresponding s-wave pion 
selfenergy diagrams appear at higher orders in density that are difficult to access through a systematic expansion of 
the selfenergy (see, e.g., Ref. [12]). In this section we will consider the renormalization of the Weinberg- Tomozawa 
interaction through vertex corrections. In Sec. V similar changes to the NLO isoscalar piece will be applied. From 
now on only symmetric nuclear matter will be considered. 

A. Tadpoles and off-shell contributions 

In the vacuum the vertex renormalizations can be partly absorbed in the coupling constant In the nuclear 
medium, these diagrams should be explicitly taken into account. Fig. 8 (1) shows a tadpole diagram that involves 
a four-pion nucleon vertex. In the free case, this term is accounted for implicitly through a renormalization of the 
lowest order Weinberg- Tomozawa term. However, in the medium the virtual pion can be polarized by exciting ph or 
Ah excitations and this leads to diagram (1') of Fig. 8. The difference between these two terms should be considered 
a genuine many body correction. 

A diagram with the same geometry but within a linear a model has been also proposed in Ref. [15]. The 47r2iV 
vertex in diagram (1) of Fig. 8 is obtained from the LO chiral Lagrangian with two baryons, 

C^^l=tTT[Br[r^,B]] (20) 
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with expanded up to four meson fields. 

^ (21) 



wliere $ is the standard SU{2) representation of the pion field, <i>ii = l/V^n", $12 = 7r+, $21 = tt , $22 = -l/\/27r". 
For the process 7r~n ^ 7r~n where the external pions have on-shcll momenta k, k' the diagrams (1) and (1') are given 

by 

which can be approximated by + fc'" ~ 2-y/s — Mi — Mj with Mi, Mj , E'i, E'j the masses and energies of the incoming 
and outgoing nucleons i and j. In Eq. (22) we have made the same s-wave projection as for the ordinary ttN — > -kN 
amplitude [31, 54]. The meson propagators for diagram (1) and (1') are given by 

00 

n 1 n [ ^ 2a.5^(u;,p,p) 

%)-p2_^2+,,' rf-(p0)2_^2+,, (23) 



where is the pion in-medium spectral function from Eq. (10). The contribution of the vertex correction can then 
be written as a correction to the kernel V ~*V + 5V of the Bethe-Salpeter Eq. (6) where 5V = V'^^ -* — V'^^\ This is 
because D{v) contains also and the vacuum diagram has to be subtracted explicitly. 

One can see from Eq. (22) that 5V has explicitly order 1//^. However, in Eq. (23) -D(i') — £'(1) is of order 1//^ 
(and higher from ph, A/i iterations in the spectral function St^) since the ph excitation p-wave pion selfencrgy is of 
order 1//^. Thus, the correction 5V is of order 1//^ and higher. 

By looking at other transitions such as 7r~p — > 7r~p or tt~p — > tt'^ti, we observe that the vertex contributions from 
Eqs. (20, 21) are of isovector nature. This means that one can absorb the vertex correction as a common factor in 
the definition of J^, as it appears in the isovector amplitude, resulting in an in-medium renormalizcd ^-^^^ for the 
isovector term, 

00 00 
-l + ^/^f-^+/rf-M-,P,p)) (24) 







with r ~ —5/12 and ~ p^ + mJ. In Eq. (24), b*i{p) and foifrcc are the density dependent isovector term and 
vacuum isovector term, respectively. Although the diagrams (1) and (1') are linearly divergent, their difference, which 
gives the medium correction, is not; thus, the p-integration in Eq. (24) is well defined. Note that casting the vertex 
correction as a correction to the coupling in Eq. (24) is just for convenience. E.g. for the ttNN p-wave coupling, 
where also appears, such a procedure does not apply. Hence, the warning here is that one should be careful not to 
talk about a universal rcnormalization of /jr. It is worth noting that for p = 0, the vacuum spectral function is given 

by 

St,{uj,P,p) 6{uj ~T]). (25) 

We observe that the integral in Eq. (24) indeed vanishes for p = 0. 

Next, we turn to another kind of contribution. In the on-shell reduction scheme of the ttN amplitude from Ref. [31] 
the on-shell and off-shell part of the ttN loop is separated and it can be shown that the off-shell part can be absorbed 
in the coupling of the TriV interaction [49]. However, in the nuclear medium, this is no longer the case and one has 
to take the off-shell part explicitly into account. In the free case, the off-shell part in the vertices of the rescattering 
diagram (2) and the crossed diagram (3) in Fig. 9 cancel the intermediate nucleon propagator in the heavy baryon 
limit, leading to a diagram with the same structure as (1) in Fig. 8. As an example we consider n^n TT~n scattering 
via a Tr~n loop as shown in diagram (2) in Fig. 9. The amplitude is then, with k ~ k' (p) the momentum of the 
external 7r~ (external neutron) and q the momentum of the tt" in the loop, 

y(^) ^j^M^/Awy I I 

J (2^)4 E{q) \ 4/2 J k'>+pO-q°- E{q) + ie - ml + ie' ^ ' 
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(2, off-shell) 



TT-k V 



\ I 



\ / TT-k' 



n 



TT-k N I ' k' 



(3, off-shell) 



/ n 



TT-k q 



r 



(4) (5) 

FIG. 9: Additional medium renormalizations at 1//^ and higher. OfF-shell parts of direct and crossed term are indicated in 
diagrams (2) and (3). Renormalization of the pion propagator in (4) and additional vertex correction with a loop in the t- 
channel, diagram (5). The shaded circles indicate resummed insertions of ph, Ah pion p-wave selfenergies in the pion propagator, 
including also A'^A^, A'A, A A short range correlations (SRC). 



Using in the heavy baryon approach — E{q) ^ and expanding the numerator as (2fc° + — fc") = 4(fc")^ + 
4fc°((7" — k'^) + ((j° — fc°)^, the on-shell part is given by the 4(fc°)^-term. For the other terms, the baryon propagator 
is canceled and Eq. (26) reads 

with V^n the usual on-shell transition TT~n — > TT~n and Gttn the 7r~n loop function. In Sec. II A medium corrections 
have been applied to the first term in Eq. (27), the on-shell one-loop rescattering. The second term is the product of 
the usual ttN on-shell amplitude times a pion tadpole and has, thus, the same structure as diagram (1) of Fig. 8. The 
remaining pion tadpole is dressed in the way it is done for diagram (1') of Fig. 8 and the vacuum tadpole is subtracted; 
the result can again be expressed in a renormalization of /^r in Eq. (24), this time with r = —3/8. The crossed term in 
TT^n TT^n scattering via one loop is displayed in Fig. 9, (3). Note that the intermediate states are in this case n^n 
and TT^p. Evaluating the off-shell parts as before, again the structure of tadpole and on-shell scattering of diagram 
(1') is obtained. Summing both off-shell parts from diagrams (2) and (3) the result can be cast in a modification of 
/tt as in Eq. (24) with r = -1-3/4. The calculation is repeated for the other coupled channels tt~p tt^p, tt~p 7r°n, 
and TT^n — s- tt'^ti and it is interesting to note that the off-shell parts of the one-loop amplitude have pure isovector 
character. This is in contrast to the on-shell one-loop amplitude with two pure isovector scatterings that results in a 
mixture of isovector and isoscalar contributions. 

In addition we have to consider structures as in Fig. 9 (4), (5) at the same order in and density. For the tadpole 
pion selfenergy in diagram (4) of Fig. 9 we consider the process tt^ti —>■ n^n with the external pions at momentum 
k and the tttt vertex given by the LO chiral Lagrangian. The vacuum tt^ selfenergy of this external pion line consists 
in charged and neutral pion loops and can be written as 

6/^ J (27r)4 ' p^ - mi 

We have written the momentum structure from the tttt vertex in a form where it becomes visible that the first and third 
term contribute to the pion wave function renormalization. These terms can be incorporated in the free pion mass. 
Continuing with the free case, the pion tadpole (4) can appear attached to an intermediate pion in the rescattering 
scheme (see Fig. 1). In this case, one of the two intermediate pion propagators of momentum k cancels the k^ — 
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structure of the second term in Eq. (28). As a consequence, a tadpole attached to the TriV- vertex results, with the 
structure of diagram (1) in Fig. 8. The medium corrections arise then from the dressing of the pion as displayed in 
diagram (1'). 

However, the pion tadpole (4) from Fig. 9 can also appear in an external pion line of the rescattcring displayed in 
Fig. 1. In this case, the first and third term of Eq. (28) contribute to the external pion wave function renormalization 
in the medium. In other words, this is a reducible diagram, because two pieces are separated by a pion propagator. 
In the search for pion sclfenergy terms we must only look for irreducible diagrams. However, the second term in the 
bracket is special because it exactly cancels the pion propagator (fc^ — toJ)~^ leading to a genuine irreducible diagram, 
that must be taken into account and is of the tadpole type of Fig. 8 (1'). 

Inserting the pion tadpole in this way in internal as well as external pion lines, the corresponding SV^^"^ from Fig. 
9 (4) is given by 

- (H) ^ "'-) Wi'I^^ - ^<«' 

which, by analogy to the terms calculated before, can be recast into a renormalization of /^r (for the purpose of the 
isovector term) given in Eq. (24) with r = 2/3. In this case the isovector character is obvious as the pion selfenergy 
is the same for all charge states of the pion. 

Note that there should be a symmetry factor of 2 as one can insert the pion tadpole also at the other pion line in 
diagram (4). However, if the pion selfenergy is inserted in an intermediate ttN loop of the rescattering series, this 
symmetry factor is not present — each intermediate pion has only one pion selfenergy insertion. Note that for the 
contribution from inserting the pion tadpole in an external pion line of the rescattering scheme of Fig. 1 there is a 
factor 1/2 to be taken into account in the wave function renormalization when considering the adiabatic introduction 
of the interaction [62]. Considering this, it is easy to see that Eq. (29) takes already correctly into account all 
multiplicity factors. 

B. Loop corrections in the t-channel 

For the vertex correction (5) in Fig. 9 we consider the process Tr^p n^p. The loop of the vertex correction is 
charged, because a neutral pion in the loop can not couple to the Weinberg- Tomozawa term. The diagram will be 
evaluated for forward scattering k = q which simplifies the calculation — this kind of approximation will be made 
several times in the following and is discussed below. Then, the vacuum amplitude for n^p —^ n^p is given by 

The pion propagators D[p) are given by Eq. (23). The term p° comes from the Weinberg- Tomozawa vertex and the 
momentum structure in the brackets is from the tttt vertex with momentum q for the external pions. By symmetric 
integration, the only non-vanishing structure is given by the combination ^{p^y^q^ . As an explicit calculation shows, 
this one-loop correction is again of isovector type. 

For the vertex correction in the nuclear medium, one of the propagators D in Eq. (30) is dressed according to 
— > with the definitions from Eq. (23). Then, a factor of two is supplied in order to account for 

the two possibilities to insert the medium dressing in either of the intermediate propagators. We have checked this 
approximation by performing the full calculation with medium dressings in both propagators which leads to negligible 
higher order corrections. 

A straight evaluation of diagram (5) along the lines used for the other diagrams gives a contribution of the same 
size as the others. However, it is easy to see that this would be a gross ovcrestimation. The contribution from diagram 
(5) can be disregarded as shown in the following. In the first place, the Weinberg- Tomozawa term, extrapolated to 
the high energies p° involved in the loop, is a gross ovcrestimation of the actual isovector ivN amplitude. In order 
to quantify this, the isovector amplitude from the unitary coupled channel model has been evaluated, with the set 
of parameters given by the right column of Tab. II. This provides a realistic amplitude up to the region of the 
A'^*(1535). Second, the (momentum dependent) ratio of the isovector amplitude of this model to the amplitude from 
the Weinberg- Tomozawa interaction has been determined. This ratio can be well approximated by a scale factor 



(31) 
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FIG. 10: Loop in the t-channel with a different p/i-insertion. 




FIG. 11: Additional vertex corrections. The two diagrams come with a relative minus sign due to the isovector character of 
the ttN s-wave interaction. 

with p the cm. ttN three-momentum, Ai = 225 MeV and A2 = 200 MeV. Including this scale factor, the medium 
contribution from diagram (5) is given by 



where we have additionally divided by the Weinberg- Tomozawa term Vwt- 

We have calculated diagram (5) in the forward scattering limit, and frirthermore assuming the external pions at rest. 
This overestimates the contribution as it is easy to see: consider diagram (5) where one external pion line corresponds 
to the pion at rest in the nuclear medium whereas the other pion line is in a rescattering loop. Then, one pion is 
at momentum (m7r,0) whereas the pion in the loop can take high values of ^^jQ. For a typical loop momentum of 
|q| = 1 GeV this momentum mismatch leads to a reduction of a factor of around two for the diagram. Together with 
the scale factor from Eq. (31), this additional reduction renders the contribution from diagram (5) small, smaller 
than the theoretical uncertainties which will be summarized in Sec. IV F. Thus, diagram (5) is neglected. 

In Fig. 10 another in-medium diagram is shown with a ph or Ah directly coupling to three pions. The SnNN- 
interaction is obtained from an expansion of the part with D, F of the LO chiral t:N Lagrangian up to three pion 
fields [64, 65]. This interaction provides terms of the form cr ■ p and cr ■ q with the three-momenta p, q of the loop 
and the external pions. respectively. As the leading term of the Weinberg- Tomozawa interaction is of the form p*^, the 
diagram vanishes by symmetric integration. 

In the following we would like to discuss another type of loop corrections in the t-channel, the diagrams of Fig. 11. 
The sum of the two diagrams involves the contribution Ud{q^ + fc") — Ud{q^ — fc"), with Ud the Lindhard function for 
only forward going bubbles. Terms involving this combination are found very small in Appendix B of [63] and we do 
not consider them. 

We do not consider selfenergy insertions in the nucleon lines. The reason is that summing over occupied states in 
Eq. (5) corresponds to a ph excitation; a local selfenergy in the particle and the hole lines cancels in the ph propagator. 
We will come back to this question in Sec. VI. 



Another renormalization of the isovector amplitude is shown in Fig. 12. The diagram exhibits two p-wave ttNN 
vertices; the same, important, short-range correlations between ph and Ah, that are included in the dressed pion 
propagator (see Eq. (11)), should also be taken into account between a p-wave vertex of the diagram and the adjoint 
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Vertex corrections from ttNN and nN A related terms 
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FIG. 12: Additional vertex correction. Dressing the pion and introducing Pauli blocking for the intermediate nucleons gives a 
density dependent correction of the isovector amplitude. 



— - — f — ■ > 

FIG. 13: The A as intermediate baryon in 7r~p n^p scattering. The tt A- vertex is in s-wave and taken from [67]. 



ph or Ah insertion in the pion propagator. The inclusion of these short range correlations (SRC) is most easily 
achieved by decomposing the pion selfenergy in a longitudinal and a transversal part Vi and Vt- The matter part of 
the diagram after subtracting the vacuum loop, divided by the tree level Weinberg- Tomozawa term Vwt, is given by 



(12) 



SRC 



Vwt 



IpI - kp) 



N 



^(/,p) 



1 



£'w(p) + ie 

2Ft(pO,p) 



-e{kF^\^\)[{Vi{p'^,Yy)'Vt{p\p) 



l-C/(pO,p)l^z(p",P) l-U{p^,p)Vt{p^v) 



-yz(/,p)-2Ft(p",p) 



(33) 



where 



^/(/,P) 



NN 



{poy 



Vt{p",p) 



firNN 



F'{\P\)9' 



(34) 



and U is the sum oi ph and Ah Lindhard functions, U = C/jv + {f-irNA/ J-wNn)'^ Ua- The monopole form factor F 
and the Migdal parameter g' have been already defined in Sec. II B. The term with 9{kF ^ p) accounts for the small 
correction from Pauli blocking of the intermediate nucleons without any modification of the pion, whereas the term 
with 9{p — kp) comes from diagrams with pion polarization through ph and Ah insertions. 

The p'^-integration is performed numerically. There is one technical complication resulting from the non-analyticity 
of the A-width in the Ah Lindhard function (step function <^{^/s — Mn — m^r); see Appendix of [59]). This leads 
to unphysical imaginary parts in 6V from the p*^-integration; the A-width is, thus, set to zero for this diagram. The 
additional short-range correlations reduce the contribution from the diagram strongly. This is in agreement with 
findings from Ref. [64] in the study of similar in-medium corrections for the isoscalar NN interaction. 

The intermediate nucleons in the diagram of Fig. 12 can also be excited. Close to threshold, even if it is off-shell, 
the A(1232) is important as we will see. The corresponding vertex correction is shown in Fig. 13. For the ttA ttA 
interaction we take the Weinberg- Tomozawa interaction of isovector type from Ref. [66] in the s-wave approximation 
of Ref. [67]. 

As in case of the corresponding diagram with nucleons discussed above, the introduction of additional SRC for the 
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(on-shell + off-shell) (on-sliell) 

FIG. 14: Triangle diagrams. The labels "off-shell" and "on-shell" refer to the Trvr-vertex. Diagram (b) is complemented with a 
diagram with the p/i-insertion in the other internal pion line. Diagram (c) is complemented with a diagram that has the loop 
on the other side of the nN vertex. 



two ttNA vertices is important; using the same projection technique as above, the result reads 



(5K 



(13) 
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TjVA 
JttNN 



Mn -p"-Ea 



1 - u{p^ p)Vi{p^, p) 1 - c/(p", p)y* (po, p; 



-Vi{p\p)-2Vt{p\v) 



(35) 



with Vt, Vi from Eq. (34), M7v(Ma) is the nucleon (A) mass, F is the monopolc form factor for the off-shell pions 
at the TriVA vertices, and /ttAta = 2.13 is the strong coupling of A to ttN . An explicit evaluation of different charge 
states shows that the corrections from the diagrams of Fig. 12 and 13 are of pure isovector nature. 

As before, the additional short-range correlations reduce the contribution from this diagram. The ratio (^VgRQ /VvvT 
can reach values up to 0.4 at p ~ po (and still 20 % at p ~ po/2). A slight modification of Eq. (35) comes from the 
fact that the ttA^A vertex, S • p, imphcitly used in Eq. (35), is defined in the A rest frame; one has to boost the pion 
momentum to this frame. This leads to a reduction of the contribution by a factor of 0.68; yet the correction is large. 

Contrary to the case of the ttA^ scattering where we shall be able to get also many body vertex corrections to the 
isoscalar part from next to leading order terms (see Sec. V), we have no control over the isoscalar ttA interaction 
which accompanies the isovector one. It is, thus, inconsistent to consider only the isovector ttA interaction as done 
here to see the effects on the ttA^ isoscalar part through rescattering. Thus, we shall not include this corrections when 
evaluating the isoscalar part coming from tiN rescattering and shall bear in mind that we have uncertainties from 
this source in the ■nA isoscalar optical potential. 

Additionally, there is a momentum mismatch if one of the external pion lines of Fig. 13 corresponds to an external 
pion at rest whereas the other pion line is inside a rescattering loop. This mismatch affects the intermediate A 
propagators; for a typical loop momentum of q = 1 GeV, the diagram is reduced by another factor of about 2.5. 
Altogether, we attribute to this source an increase of hi from 10 to 20 % at p = po: accepting this band as theoretical 
uncertainty. 



D. Triangle diagrams 

There is another family of diagrams displayed in Fig. 14. As indicated in the figure, it is enough to calculate the 
diagram on the right hand side with the tttt vertex taken at its on-shell value in the sense that, whenever p^ appears 
in the tttt amplitude, it has to be replaced by m^. This is equivalent to calculating the same diagram on-shell plus 
off-shell, plus the set of other diagrams displayed on the left hand side of Fig. 14. This has has been shown in [64, 65] 
and has been verified again in the present study for the limit of zero momentum exchange which is also the limit 
taken here for all vertex corrections. Such cancellations were first shown by [68, 69] in the problem of tttt scattering 
in the nuclear medium. 

For the calculation, we consider first the reaction n^p —>■ ■n^p. For this configuration of external particles we 
can have charged or neutral pions for the loop lines. Summing both possibilities and inserting a factor of two from 
inserting the medium correction in either internal pion line, the medium amplitude takes the form 
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FIG. 15: Additional set of triangle diagrams with A(1232). The figure caption from Fig. 14 also applies here. 



where the momentum of the external (internal) pions is k (p), respectively. Again, we take the limit of forward 
scattering and, moreover, that the external pions and nucleons are at rest. In Eq. (36), E{p) = a/ + |pp is the 
nuclcon energy. There is also a form factor -F(|p|) for the off-shell pions in the ttNN vertex and the factor /E{p) 
from the non-relativistic reduction of the nucleon propagator. 

A straightforward calculation reveals that the term 3p^k^ is of isovector nature, whereas the contribution with 
3/4to^ is isoscalar. In the heavy baryon approximation we can neglect Mjv — E{p) in the baryon propagator and in 
this limit the isoscalar term cancels due to symmetric integration. We are left with a purely isovector contribution in 
which the p'^ from the numerator cancels the baryon propagator in the heavy baryon limit; the correction is given by 



1 , f ^ S,(u>,P.p) 



4r] 



2 



(37) 

The term with 9{kF ^p) accounts for the small medium correction from Pauli blocking of the intermediate nucleon but 
without any modification of the pion, whereas the term with 9{p — kp) contains all diagrams with pion polarization. 

In Fig. 15 we plot the analogous diagrams of Fig. 14 but with a A intermediate state, instead of a nucleon. The 
same type of off-shell cancellation found for the diagrams in Fig. 14 holds also here as has been shown in [64, 65] for a 
similar configuration. This means, with p (fc) being the loop momentum (external momentum): p^ mj, fc^ mj, 
pk — > (the term pk is not affected by the off-shell cancellation). In the last substitution, the integration over the 
spatial part pk vanishes. 

Including a factor of two from inserting the medium correction in cither pion line, the correction for n^p ir^p is 

oo 

OO 

ml{MN -Ea- 2r,) - 2if k\ [ ^ S^{io,p,p) 2ml (Mn - E^ ~ lu ~ j^) - ^lutj fc° ^^^^ 



^fifni J 2tt^ t] 







2if{MN - Ea-v)^ J uj + rj [Mn - EA-ri)(MN - Ea-oj) 



Here, Ea is the A(1232) energy. As an explicit calculation shows, the term with 2m^ inside the w-integral is of 
isoscalar nature whereas the term with k'^ is of isovector nature. This means that the term with 2mJ is the same for 
all channels of our coupled channel approach whereas one has to multiply SV^^^^' '^^ by —1, — \/2, and for ir'n 7r~n, 
Tr~p —> Tr'^n, and ir'^n Tr^n, respectively. As nucleon and A mass arc non-degenerate, the isoscalar part docs not 
cancel as it had been the case for the diagrams of Fig. 14. 

In Eq. (33) we have already seen an example for additional SRC: Between the nucleon emitting a pion and the 
ph or Ah medium insertions in the pion propagator, there are also SRC. If this is the case for both ends of the pion 
propagator, the substitution of the pion propagator D^j is given by the projection technique employed in Eq. (33). 
In the triangle diagrams, there is only one side of the pion line affected and the additional SRC can be cast in a 
substitution of the in-medium pion propagator -D(i') from Eq. (23) according to 

D^,,)^D^y) X (39) 

as an explicit calculation shows. This does not affect the off-shell cancellation behavior discussed before. Remember 
that we always use the same Migdal parameter for NN and NA SRC. Although the contributions from Eqs. (37,38) 
are sizable, they are suppressed by a factor of five by the additional SRC from Eq. (39). In the final numerical results 
they introduce a small correction. 
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FIG. 16: Nucleon tadpole in the pion propagator with NLO chiral -kN interaction. 

E. Isovector correction from the NLO ttJV interaction 

The Weinberg- Tomozawa term is also renormalized by higher orders in the isoscalar t^N interaction. A correction 
of this type comes from the nucleon tadpole in the pion propagator shown in Fig. 16. The ixN interaction of the 
tadpole is from the s-wave isoscalar interaction from the NLO chiral Lagrangian whereas the other -kN interaction is 
given by the Weinberg- Tomozawa term. The nucleon tadpole with isovector interaction vanishes in symmetric nuclear 
matter. 

The isoscalar -kN interaction in the NLO order is given by [46] 



J 71 J TT J 7T 

. - ^ u^ikf ^ ml) - ?^i±^ (fc^ - ml) = + t%- (40) 

Here the interaction has been separated into on-shell and off-shell part [70] (the latter term with (fc^ — ml)). For the 
nucleon tadpole in Fig. 16 and considering first the off-shell part, the selfenergy is given by 11 = — (2c2 + 2c3)(fc^ — 
ml)p/ f^. The entire diagram is then given by 

= - ml)D{k)p (41) 

Jtt 

where tT^M^-^N is the isovector interaction from the Weinberg- Tomozawa term and D{k) the intermediate pion prop- 
agator which cancels the term (fc^ — ml) from the isoscalar vertex. This means a vertex renormalization by a similar 
mechanism as we have already seen for diagram (4) in Fig. 9 and 

5V^'''^ 2c2 + 2c3 

1^'- ^''^ 

For the numerical evaluation we use the values of the c-cocfficients of the fit 2'^ from Ref. [52], 

ci = -0.35 ±0.1 GeV"\ ca = -1.49 ± 0.67 GeV"\ = 0.93 ± 0.87 GeV"\ (43) 

It would be more consistent to use the values of the present fit in Tab. II instead. However, in the present model we 
have only access to C2 and the combination 2ci — C3. In Eq. (42) the c-coefficients are combined in a different way, 
and we have to resort to the values of [52]. In any case the values from Eq. (43) are compatible within errors with 
ours from Tab. II (see also Ref. [31]). For an estimate of the theoretical error, we can also use the c- values from fit 2* 
instead of 2^ [52]. This induces a theoretical error of the order of 20 % for the contribution which by itself is smaller 
than other diagrams. 

As for the on-shell part of the interaction in Eq. (40) we notice that the intermediate pion propagator, between the 
nucleon tadpole and the Weinberg- Tomozawa vertex, does not cancel. This means that the on-shell nucleon tadpole 
contributes to the pion selfenergy and not to the vertex renormalization and hence, if it is an external pion line, it 
will be automatically considered when solving the Klein Gordon equation for pions with the proper selfenergy and 
must not be included as a genuine new contribution. When this part of the tadpole occurs in internal lines, compared 
to the p-wave pion selfenergy, the s-wave selfenergy is small and can be neglected as we have also seen in the self 
consistent calculation in Sec. HI A. 

F. Results for the isovector renormahzation 

For all corrections evaluated in this section IV, the vertex corrections can be recast as a correction to the isovector 
interaction strength h\{p) or, in other words, an in-medium change of /t^. Note that we refer to the /^r that appears 
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FIG. 17: In-medium isovector bl{p) compared to the vacuum isovector term feifrco- The gray band from Suzuki et al. [26] is 
from a phenomenological fit, as well as the point from Nieves et al. [22]. Also shown are chiral calculations from Meifiner et al. 
[25] and Weise et aL (Friedman et al.) [13, 71]. 



in the Weinberg- Tomozawa term of Eq. (3); wc do not claim a universal change of f-^ in the nuclear medium (see also 
the caveat following Eq. (24)). For example, Eq. (24) gives the renormalization of feivac/&i(p) from the diagrams (1) 
to (4) of Figs. 8,9 with an overall value of r = 1. Including these diagram as well as all other isovector corrections 
found, the in-medium change of bi is plotted in Fig. 17. 

The result in Fig. 17 is given at A = 0.9 GeV for the monopole form factor that appears in the ph and Ah pion 
selfencrgies of the vertex corrections (see Eqs. (11,34,39)). The dependence on A is moderate and at A = 1 GeV 
^i(p)/^i,frcc is increased by another 10 %. With the decrease of in the medium the isovector ttN interaction 
effectively increases, in quantitative agreement with a recent analysis of deeply bound pionic atoms [26] and the 
phenomenological fit from [22] and the analysis of ttA scattering from [71]. There is also qualitative agreement with 
other theoretical works [13, 25] which are justified in a different way, without the thorough many body study done 
here. 

In Fig. 17 we give a band of values for our results of bi, including the uncertainties discussed above from the 
diagram of Fig. 13, plus a 20 % extra uncertainty from the dependence on the form factor. We can see that the band 
overlaps with the experimental band of [26]. 

After studying vertex corrections of isovector type in this section, next we turn to study vertex corrections of 
isoscalar type and their effects in the isoscalar optical potential. 



V. RENORMALIZATION OF THE NLO ISOSCALAR TERM IN TriV SCATTERING 

The model from Rcf. [31] for the ttN interaction in the vacuum has two sources for isoscalar contributions: one 
is the NLO, point-like, interaction from Eq. (4) and the other one comes from the isovector term in the rescattering 
of the pion generated in the Bethe-Salpeter equation. In fact, the latter is quite large, bg = 442 • 10^'*7ti^^ (see 
Tab. VII from [31]). This large contribution is partly canceled by the NLO contact term from Eq. (4) that is 
be ~ —336 • 10^'' ^^^^^7 leading to a final value of bo = —28 • 10^* "i^^. 

For the application of the model in nuclear matter this partial cancellation has consequences. Renormalizing the 
isovector strength changes the in-mcdium isoscalar term through the rescattering piece. Then, the sum of this term 
and the point-like NLO interaction will not show the partial cancellation of the vacuum case any more. It is therefore 
important to treat the NLO isoscalar term on the same footing as the isovector renormalization. 

It is easy to see what the effects of the increase of bi will be in the isoscalar part of the potential. Indeed, if one 
cared only about Pauli blocking corrections in the intermediate nucleons, the effect would be given by Eq. (16) with 
an increased 6i, thus leading to an increased isoscalar repulsion. However, the TriV rescattering term with no Pauli 
blocking, which is larger, has opposite sign (see sign of bg). Hence, the net effect of increasing bi in the medium, 
including Pauli blocking corrections, is a net attraction, with opposite results to a naive implementation of the bi 
changes in the Ericson-Ericson formula of Eq. (16). In some analyses of data [71] the needed repulsion is obtained 
by using the Ericson-Ericson Pauli blocking correction, Eq. (16), with the increased bi. This is one empirical way to 
implement repulsion, but from the theoretical point of view one should evaluate the full ttN t-matrix in the medium, 
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FIG. 18: Pion-pion interaction with a nucleon tadpole Irom the NLO vrA'^ interaction (left). Closing one pion line produces a 
pion selfenergy of the tp type (center). This t is a vrA*' vertex correction (right), with the same geometry as diagram (1) in Fig. 
8. 

not just the Pauli correction to it, and this means that the problem of the missing repulsion becomes more acute. 

The diagrams from the last section will be the guideline for the renormalization of the NLO isoscalar. We do not 
redraw these diagrams, but the t:N — s- ttN contact interaction is now given by the NLO isoscalar term instead of the 
LO Weinberg- Tomozawa isovector term. 

A. Tadpole and off-shell contributions 

We start with the pion tadpole (1) from Fig. 8. The NLO isoscalar Lagrangian has to be expanded to four pion 
lines in order to provide the 4:1t2N vertex needed in this diagram. As shown in the following, to this end we can utilize 
the in- medium Lagrangian derived in Ref. [72-75] (see also [70, 76]) by taking the mean- field approximation for the 
nucleon field. The terms with a medium correction p of the nuclear density read 

(£) = i p (C3 Tv[dUdU^ + C2 Tv[d„Ud„U^ + Tt[U^x + X^U]) (44) 

by keeping only the isoscalar terms which are parametrized in terms of ci, C2, C3. Expanding this term up to four 
external pion lines leads to a tttt vertex with a nucleon tadpole as displayed in Fig. 18 to the left. Contracting two 
of the pion fields leads to a diagram that appears as a pion selfenergy with a pion tadpole and a nucleon tadpole as 
displayed in the center of Fig. 18. For a tt~ this selfenergy is given by 

n(^^) = ^*/ (04^W (2c3(fc^+p^) + 2c2((fcY + (pY)-^crm^) (45) 

with the pion propagator D{p) and the momenta as assigned in Fig. 18, center. The Lagrangian of Eq. (44), as 
well as its related one of Eq. (4), are meant to be used at very low pion energies and their extrapolation to high 
momenta is not justified. Actually, in the study of ttN of Ref. [31], where an extrapolation to energies about 400 
MeV above ttN threshold was done, the exponential damping factor of Eq. (4) was demanded by a fit to the data. 
Hence, we make here the sensible choice of setting = (p°)^ = in Eq. (45). While this certainly introduces some 
uncertainties, these are still small compared with larger sources of uncertainties that we shall discuss below. 

Note the appearance of the nuclear density p in Eq. (45): the selfenergy is of the type tp with a matrix element 
t that can be extracted by opening the nucleon line of the nucleon tadpole, meaning the division of Eq. (45) by p. 
This is displayed to the right in Fig. 18. The resulting diagram is a irN vertex correction with the same geometry as 
diagram (1) in Fig. 8 but using the NLO isoscalar interaction for the 4:Tt2N vertex. 

Next, we proceed like in the diagrams of Fig. 8 including ph, Ah, and short range correlations in the pion propagator 
of the third diagram of Fig. 18 and subtracting the free part of this ttN matrix. This gives a genuine many body 
correction, in the line of Eq. (24) but of isoscalar character. As a result, the isoscalar vertex correction for the coupled 
channels i , j reads 

00 00 

<'' = l|/lf Jd^[sA-,P,p)-'-^) [m^(8c3+4c2-5cr) + C2(2fcO)^]. 


(46) 

The in- medium correction from Eq. (46) has to be added to the kernel of the Bethe-Salpeter equation (6). The term 
with corresponds to the dressed pion propagator as in diagram (1') of Fig. 8 and the term with S{rj — lu) to the 
vacuum diagram (1) which is subtracted. 
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For the diagrams (2) and (3) from Fig. 9, one or two of the irN vertices can be given by the NLO isoscalar 
interaction: the Bethe-Salpeter equation (6) iterates the kernel and allows for any combination of isoscalar (see Eq. 
(4)) and isovector vertices (see Eq. (3)) in the rescattering series. For two iterated Weinberg- Tomozawa vertices we 
have already determined the off-shell contributions from the direct and crossed term; for iterated isoscalar interactions 
we should do in principle the same. However, the contributions are smaller and we neglect them. This can be seen 
as following: the strength of the isoscalar interaction is be = —336 • 10^'*to^^ [31] which is around one third of the 
isovector strength. A combination of two isoscalar vertices in TriV rescattering would, thus, approximately lead to 
an off-shell contribution nine times smaller than the off-shell effect from the combination of two isovector vertices 
studied before, and we can safely neglect it. A combination of one isoscalar vertex and one isovector vertex results 
in an overall isovector interaction and is of no interest in the present case where only symmetric nuclear matter is 
considered. 

The diagram (4) from Fig. 9 renormalizes the NLO isoscalar interaction in the same way as it affected the isovector 
interaction studied before. This means SvI^Iq/Vnlo is given by the right hand side of Eq. (24), with r = 2/3 as 
before, where Vnlo is given by the second term in Eq. (4). 



B. Loop corrections in the f-channel 



In diagram (5) of Fig. 9 the Weinberg- Tomozawa term can be replaced by the NLO isoscalar interaction. We will 
consider on-shell and off-shell part of this interaction, given in Eq. (40), separately. 



On-shell part of the ttN — > ttN vertex 



We consider the process tt p n p and a charged pion running in the loop. Then, the contribution from the 
on-shcll part of the NLO Lagrangian reads (sec first term of the right hand side of Eq. (40)) 



-zT, 



(5), on _ 



NLO 



d'^p 

— — - {-iUn) iD(p) iD{p) {-iU^) 
(27r)^ 

1 f _d'^p Aciml ~ 2czml - 2c2Ti^ 
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J TT 



(p2 



■ ie)^ 



[p" + &pq + q^~ 2ml) 



(47) 



with 



ml 



p ^ and D being the pion propagator as before. Wc have taken here already the limit of forward 



scattering. The last term comes from the tttt vertex with q the momentum of the external pions. As we did before. 



we take the on-shell value for the external pions q 



m' 



and can substitute the last term in Eq. (47) by {p — m^). 



taking into account that the mixed term Qpq vanishes due to symmetric integration. As in Sec. IV B we dress only 
one of the propagators in order to stay in line with the other corrections evaluated (then, a multiplicity factor of 2 
appears according to the two possibilities of inserting the in-medium correction in either pion propagator of the loop). 



in Eq. (47) for the medium part and D 



^(1) 



for 



The propagators in Eq. (47) are then given by D L'(i) Z?(i/ 
the vacuum part with £*(!), from Eq. (23). 

The NLO ttN amplitude in Eq. (47) has a term 2c277^ which introduces additional powers of p in the integration. 
Although the integral is still convergent, the value of the C2 coefficient is only valid for small momenta, where it has 
been determined in fits to low energy irN scattering data (see Eq. (43)). Once more, in analogy to what was done 
following Eq. (45), we replace 2c2ry^ ~> 2c2mJ in Eq. (47), taking thus the threshold value of the NLO isoscalar 
interaction. 

Furthermore, there is another off-shell cancellation, the one of the tttt vertex. This is due to the additional diagram 
shown in Fig. 10. In a similar way as in [64, 65] (sec also the diagrams in Fig. 14), the diagram cancels the off-shell 
part of the tttt interaction appearing in Eq. (47). 

For the process n^p 7i'~P, the pion in the loop can also be a 7r°. It is easy to see that the overall correction, 
including all the intermediate pions with different charge, is of isoscalar type by explicitly calculating other ttN 
channels. Taking all this into account, integrating the p°-component, and subtracting the vacuum part from the 
medium part results in 
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(48) 



This is a tiny correction to the isoscalar rcnormalization which is neglected in the final numerical results. 
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2. Off-shell part of the ttN — > ttN vertex 



The off-sheU part of the NLO isoscalar interaction is renormahzed in a similar way as before. Taking the term with 
^ — of Eq. (40), the vacuum amplitude for ir^p n^p is in this case given by 



(5),off 
NLO 



2c2 + 2c3 1 



/2 3/2 y (2^)4 (p2_^2 + 



(p^ - ml) 



(p2 + 6pq + - 2ml) + ip^ + 9^ 



(49) 



where the external pions are again at momentum q and (p^ 



is the off-shell part of the NLO isoscalar vertex 



which cancels one of the propagators. In the square brackets, the contributions from having a charged pion or a 
neutral one are denoted separately. It is easy to see that the overall contribution is again of isoscalar nature. 

In the term (jp + 6pq + q^ — 2ml) from the tttt vertex, the mixed product pq vanishes due to symmetric integration. 
The ^''-integration is straightforward for both vacuum and medium loop. The resulting medium correction reads 
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(5),off 
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dp p^ 
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+ dio Sn{Lo,p, p) u 



(50) 



The factor (oj^ 



p^ 



il j 4) comes from on- and off-shell part of the tttt- vertex. As before, the off-shell part cancels 



with the diagram from Fig. 10. Keeping only the on-shell part of the tttt interaction, the final result can be easily 
obtained and is given in Eq. (50) with the replacement 



2 2^21 3 2 

w -p - -m^j ^ -m^ 



(51) 



C. Further renormalizations of the isoscalar nN interaction 



Next, the diagram from Fig. 12 is considered which is now given by two ttNN vertices and one Tr~N Tr~ N 
transition from the NLO isoscalar interaction. The only change with respect to the previous result from Eq. (33) is 
a change in isospin factors. Indeed, consider the n^p n^p amplitude and the loops with ir^p or 7r+n (see Fig. 12), 
implying the ttN — > ttN vertex to be Tr~p — > n^p or n^n — > Tr~n, respectively. The isovector TriV vertex has opposite 
sign in these two cases while it is the same one for the isoscalar amplitude. Hence, the contributions of the two loops 
add in the case of the isoscalar correction while they get subtracted in the case of the isovector correction. Thus we 

obtain (JV^^lq/^nlo = — 3 JVg^^V^T with SV^^^q /Vwt given by Eq. (33). 

Considering the nuclcon tadpole in Fig. 16 it is clear that this correction renormalizes in the same way the isoscalar 
and the isovector interaction because the tadpole factorizes with the t:N amplitude. Thus, the renormalization of the 
isoscalar amplitude (^Vj^lq/V^nlo is again given by the right hand side of Eq. (42). 

When renormalizing the Weinberg- Tomozawa interaction, we have also considered the diagrams in Fig. 11. They 
have been found small as argued at the end of Sec. IV A due to the occurrence of the difference of Lindhard functions 
U{q^ + fc°) — tj{q^ — fc°). The minus sign was a consequence of the isovector nature of the Weinberg- Tomozawa vertex. 
However, if one replaces the Weinberg- Tomozawa vertices at the bottom of the diagrams in Fig. 11 with isoscalar 
ones, one obtains the combination U{q^ + kP) + tj{q^ — fc°). This might result in a significant contribution. However, 
as the following argument shows, we should not consider this contribution as it would mean double counting: We 
consider the diagrams in Fig. 11 with the nucleon line closed. This corresponds to a contribution to the s-wave 
pion selfenergy. However, the resulting selfenergy, let it be H(ii), is already generated in a different piece: Imagine 
the one- loop rescattering of two isovector interactions. Insert now a nucleon tadpole in the intermediate pion line as 
displayed in Fig. 16. Joining the external nucleon lines in this rescattering diagram generates again the pion selfenergy 
H(]^x)- For similar reasons of double counting we also discard vertex corrections that occur when one replaces the 
vertices at the bottom of Fig. 11 with the NLO isoscalar interaction and, additionally, also the other -kN vertices. 

Finally, let us remember there is a small vertex correction to the isoscalar interaction from the triangle diagrams 
with intermediate A of Fig. 15 which we have already discussed and evaluated in Sec. IV D, Eqs. (38,39). 

The sum of the isoscalar corrections calculated in this section results in an increase of the NLO interaction as shown 
in Fig. 19. On the right hand side, the 56^'"'-*, obtained from the various corrections 6Vnlo of this section, is shown, 
divided by 6o from the vacuum NLO isoscalar term given by the second term of Eq. (4). The ratio is of similar size 
as the isovector case shown in Fig. 17. On the left hand side of Fig. 19 we show again the NLO isoscalar 6o term 
for vacuum and medium. Additionally, the vacuum isoscalar term from the rescattering of isovector interactions is 
plotted, indicated with "&o rescattering (vacuum)". This illustrates how these two sources of isoscalar contribution 
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FIG. 19: In-medium isoscalar interaction. Left: different sources of isoscalar contributions are shown (see text). Right: The 
fraction 6bo^'^ /bo^'^ , actually calculated here, is plotted as a function of kp. 



almost cancel in the vacuum. When the full ttN vertices arc used in the rescattering, the small "final bo (vacuum)" 
appears, as required by the vacuum data [31]. Fig. 19 also shows, that the isoscalar in- medium vertex corrections, 
evaluated in this section, provide a large source of repulsion (negative 6o), even larger than the Ericson-Ericson 
correction from Eq. (16), which is also shown in the figure. In Sec. VII A we will test the final results of the medium 
calculation for changes of the vacuum model. In particular, we will perform refits requiring a smaller NLO isoscalar 
vacuum term, which then also reduces the size of isoscalar contribution from rescattering (see Fig. 19). However, as 
will be seen in Sec. VII A, the final results are stable under these variations. 



VI. METHODOLOGY OF THE EXPANSION 



The diagrams introduced follow the standard approach of field theory in which selfenergy corrections as well as 
vertex corrections are introduced in a perturbative expansion. Here, we want to discuss the methodology behind the 
expansion done. Let us proceed step by step. The first step consisted in developing an accurate model for the free 
ttN scattering amplitude in Sec. II. There, we have followed the standard procedure of the chiral unitary approach. 
The Bethe-Salpeter (BS) equation is used to ensure unitarity with a kernel derived from chiral Lagrangians in which 
the lowest order and next-to-lowest Lagrangians are considered. An expansion in powers of momenta of the pion is 
done in the kernel of the BS equation, and the BS equation resums the higher order terms, like in the Schrodinger 
equation starting with a potential, ensuring the unitarity of the amplitude. 

When this is done, the many body expansion comes into action. Attempts of a systematic many body expansion 
along standard lines of chiral perturbation theory have been made in [25], but one automatically faces problems 
ultimately tied to the fact that the many body corrections come from ph excitations which have small energy (starting 
from zero in our Fermi sea). This is in contrast with chiral perturbation theory where the success, for instance in tttt 
scattering, is tied to the large energy gap between the tt and the next meson mass excitation, the one of the p. Even 
in TT-deuteron scattering many options of expansion parameters are suggested [41] with quite different results among 
them. 

Our guiding line for the many body expansion follows the traditional many body approach in nuclear physics where 
one relies upon an expansion in powers of the density. The density is assumed to be a small quantity and this is based 
in phenomenological facts. Let us illustrate this in our case of 7r-nucleus scattering. Let us take Fig. 4 where the 
selfenergy of the pion is obtained by summing the irN T-matrix over the occupied states of the Fermi sea. Technically, 
we have introduced a hole line in the diagram which induces a power of p upon integration over the occupied states. 
Let us now cut the diagram by a vertical line, cutting the pion and the hole line. Restricting ourselves to one meson- 
baryon loop for simplicity, we have the result of Fig. 20. The whole sum of loops would replace the itttNN vertex 
by the full ttN T-matrix. Placing on-shell the intermediate states, cut by the vertical line in Fig. 20, as it would 
be implemented using Cutkowsky rules [77], one accounts for tt transitions to a tt -I- (ph) channel. In other words, 
the diagram accounts for ttN — > TriV, which is the quasielastic scattering. Let us now dress the intermediate pion 
by letting it excite a ph, like in Fig. 21, and let us now cut the two ph excitations by a vertical line placing them 
on-shell. The process that this cut accounts for is what is called genuine two-body absorption. One has ttNN NN 
and the pion disappears. Experimentally one has, for energies not close to threshold, bigger quasielastic cross sections 
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FIG. 21: Two-nucleon pion absorption. 



than absorption ones. Close to threshold, Pauli blocking makes the quasielastic scattering small, which means that 
the imaginary part of the diagram of Fig. 20 becomes small, but this is not the case for the real part which is not 
subjected to the Pauli exclusion principle [78] and which has the size of the imaginary part when Pauli blocking 
does not restrict it. Much work has been done in the past on three body absorption, both experimentally [79] and 
theoretically [80], obtaining smaller rates than for two-body absorption, particularly at energies close to threshold. 
Thus, an expansion in the number of ph lines, or in general terms, a hole line expansion, is the guiding principle for 
our many body approach, the justification for it coming from phenomenology. 

This said, let us see what a strict bookkeeping along the lines described above would give us when we introduce 
an extra ph in the basic diagrams. Start from the loop diagram of Fig. 22(a) and introduce a new ph. Fig. 22(b) 
is accounted for by our diagram of Fig. 3. Fig. 3(c) plus 3(g), external line renormalizations, introduce a selfenergy 
in the nucleon line when folded and converted into a hole line like in Fig. 23. This has to be considered together 
with diagram (e) in the sense that both particle and hole have to be renormalized simultaneously to account for the 
selfenergy in the medium. Yet, this selfenergy, ignoring non-localities, as done in our approach, gets cancelled in the 
ph propagator [Et^ + Eh — Ep]^^. The diagrams (d) and (f) have already been accounted for in the diagram of Fig. 
12. Diagrams (h), (i), (j), (k) have already been considered before in Fig. 14. Similarly, diagrams (m) and (n) are 
taken into account by the vertex correction from Fig. 8. We are left as novel diagram with the one of Fig. 3(1) which 
involves the ttttN intermediate states that were considered in [54] and found to be relevant for / = 3/2 at energies 
around s^/^ = 1.5 GeV, very small for / = 1/2 at these energies, and negligible at nN threshold that we study here, 
as we noted at the end of Sec. II A. We use this fact here to neglect the medium modifications to this negligible term 
in free space. 

The discussion presented before clarifies further the procedure followed in former sections, where following a tradi- 
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FIG. 22: Systematical ph insertions in the basic rescattering loop. 
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FIG. 23: Hole self energy. 



tional approach in field theory one studies selfenergy and vertex corrections. 

In addition to this, our approach introduced further ph excitations following again traditional lines in field theory 
and many body. For the selfenergy corrections, see Fig. 3, the ph excitation is induced as a selfenergy of the pion, 
and thus, iterated ph excitations in the Dyson sense are automatically accounted for. This is done technically with 
no extra effort and produces the coupled branches (tt and ph, or A/i) of the pion in the medium [81]. 

The self consistency that we have implemented, introducing the full calculated tt selfenergy into the intermediate 
TT of Fig. 3 also accounts for the s-wave self energy of higher order and is always desirable in many body calculations, 
although its effects here are moderate. This is unlike the case found in the i^~-nucleus interaction [55] where the 
presence of a resonance close below threshold renders this procedure rather important. 

The discussion in this section has clarified further the principles followed in the many body expansion, which are 
ultimately based on phenomenological facts and not in a formal expansion based on one single parameter, as one has 
in some theories for elementary particle interactions. 

VII. NUMERICAL RESULTS 

In the last sections IV and V, vertex corrections for both the isovector and the isoscalar interaction have been 
evaluated. Together with the in- medium ttN loops Gt^n, shown in Eq. (8), we obtain a new irN — > ttN transition 
T in Eq. (6). Integrating over the nucleons of the Fermi seas according to Eq. (5), the s-wave pion selfenergy is 
evaluated. The results are for symmetric nuclear matter as the calculations from Sees. IV and V are performed in 
this limit. 

The contributions to the s-wave pion self energy can be ordered in powers of the density. In the following, we 
discuss the contributions up to order p^p^^^ for the external pion s-wave selfenergy IVs- In Sec. VII A, we compare 
this to the result up to all orders in p. In the ■nN rescattering loops, the pion p-wave polarization accounts for one 
power of p whereas the closing of the external nucleon line, as indicated in Fig. 4, corresponds to another power of 
p. The Pauli blocking of the intermediate nucleon in the ttN rescattering generates an extra power of p^/^ as can 
be seen in Eq. (16). Of course, higher powers are also contained: from the multiple rescattering generated by the 
BSE equation (6) on one hand, and the resummation of ph, A/i pion selfenergies from Eq. (11) on the other hand; 
however, these higher order corrections are small. 

Another contribution to Us comes from the vertex diagrams from Sees. IV, V, which are of order p through the 
p-wave pion polarization. Closing the nucleon line (order p) for these diagrams, corrections of order p^ are obtained. 
The vertex diagrams induce also corrections at (P'p^l^: those consist of the Ericson-Ericson rescattering piece from 
Fig. 3 without pion polarization of the rescattered pion, but with vertex corrections from Sees. IV and V for exactly 
one of the Weinberg- Tomozawa "kN interactions. Some examples can be seen in Fig. 22. These contributions are 
most easily included by multiplying the rescattering term by (fcx(p)/&i.fioc — 1), with fci(/3)/6i.frco from Fig. 17. 

Summing all contributions, the pion s-wave self energy up to order (P' p^l^ is plotted in Fig. 24 with the black 
solid lines. For the pion three-momentum, we have taken a typical value of |p| = 50 MeV although Ilg depends only 
weakly on p. The gray band shows the area of the experimental fits to pionic atoms from Refs. [27, 28] (see also 
Fig. 6) whereas the dark band represents the phenomenological fit from [22]. The present result for the external pion 
energy fc" = ttItt stays some 30 % below the phenomenological fit. Note that at the order p^ p^l^ considered here, the 
imaginary part Imllg is the same as in Fig. 7. 

In Refs. [82], [14], [17] it has been claimed, that a possible way of understanding the repulsion in pionic atoms 
comes from the energy dependence of the pion self energy; a consistent treatment of the Coulomb potential in the 
Klein-Gordon equation requires that the argument of the pion self energy is 115(11; — V^) [14] rather than 115(0; = ra-,^). 
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FIG. 24: The s-wave pion selfenergy in nuclear matter up to order p^p^^^, for three pion energies (fc" — m,r, = m,r -I- 
lOlVIeV, = m.^ + 20 MeV). The gray band shows the area of the phenomenological fits from Refs. [27, 28] and the darlt 
band the fit from [22]. For Imlls all phenomenological fits lie in the gray shaded area. 
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FIG. 25: s-wave pion selfenergy to all orders in p. Different pion energies (fc" — m.^, fc" = + 10 MeV, fc" = + 20 MeV). 
Phenomenological fits as in Fig. 24. 



Furthermore, the small isoscalar ttN potential at threshold rises rapidly with increasing energy and, thus, large effects 
from the energy dependence of lis can be expected. 

Ill order to see this effect in the present calculation we have plotted Ilg also for p'^ = m„ + 10 IVIeV and p'^ ~ m^r + 20 
MeV (dashed and dotted line, respectively). The Coulomb potential for a nucleus with A = 100, Z ~ 50 can reach 
Vc ~ 16 JVIeV at an effective density of p = po/2, and even more for heavier nuclei. As Fig. 24 shows, the energy 
dependence leads indeed to an extra repulsion which agrees well with the phenomenological fits. However, theoretical 
uncertainties are larger than thought as will be discussed in Sees. VII A, VII B. 



A. Theoretical uncertainties 



Certain higher order corrections in density play an important role but at densities higher than felt by pionic 
atoms. E.g., at order p^p^^^, diagrams appear where both isovector vertices of the large rescattering piece contain 
the corrections from Sec. IV. These corrections reduce Re Us from Fig. 24. Dressing all vertices with the corrections 
found and using the ttN loop function with Pauli blocking and pion polarization, i.e. including all corrections found, 
to all orders, gives the result indicated in Fig. 25. Changes with respect to Fig. 24 are mainly due to the fact, that 
the vertex corrections can occur quadratically and higher, while before, at order p^p^^^, only one vertex correction 
enters the rescattering series. 
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The imaginary part on the right hand side in Fig. 25 is more negative and closer to the values from phenomenological 
fits: the imaginary part comes from the ph insertions in the ttN loop function. The imaginary part of this rescattering 
loop is enhanced by the larger strength in both ttN vertices due to the vertex corrections. 

The decrease of the real part can be understood as following: in the vacuum model the rescattering piece introduces 
an attraction which is compensated by a repulsion be from the NLO isoscalar interaction at tree level [31] . The isovector 
interaction in the medium is increased from vertex corrections as we have seen in Fig. 17; this leads to an increase of 
the attraction from the rescattering piece. The Pauli blocking of the intermediate nucleon, namely the Ericson-Ericson 
effect from Eq. (16), which is repulsive, can not fully compensate this effect; as a result, the net repulsion is smaller 
than without vertex corrections. However, for arguments of external pion energies of k'-' = 771,^ + 10 MeV, fc*^ = m^r + 20 
MeV due to Coulomb shift, the results are in the region of the required repulsion at kp ~ 210 MeV {p ~ /Oo/2). It 
should be noted that the net attraction of the the increased bi in the medium, through rescattering, is a rather large 
effect, which more than compensates for the enhanced isoscalar amplitude in the medium found in Sec. V and shown 
in Fig. 19. 

We have tested the stability of our results. The corrections discussed stemming from rescattering are tied to our 
elementary vacuum model of Ref. [31]. For instance, the size of the isoscalar contribution be from the NLO Lagrangian 
is correlated with the subtraction constant of the nN vacuum loop, OttTv, in the fit. We have performed a refit of the 
vacuum amplitude requiring a smaller value of be- Close to threshold, a sufficiently good fit can be obtained with 
be = —46 • 10~^m~^ which is around ten times smaller than be from the fit of Tab. II. It is interesting to note that 
the resulting parameter values 2ci — C3 = —1.43 GeV""'^, C2 = —1.54 GeV~^ are quite close to the values from [52] of 
2ci — C3 = —1.63 ± 0.9 GeV~^ and C2 = —1.49 ± 0.67 GeV~^. However, even with this drastic change of the vacuum 
model. Re H5 hardly changes and the results are stable in this respect. 

Further theoretical uncertainties come from the regularization scale A that appears in the monopole form factors 
of the pion p-wave polarization. The result depends on A; a smaller value than the one used of A = 0.9 GeV would 
provide slightly larger repulsion as we noted before. Nevertheless, the good agreement with the phenomenological 
analysis on the 61 renormalization from [26], which has been noted in Sec. IV F, provides support for this value of A. 

We have also treated the pion self energy selfconsistent as in Sec. Ill A, including the energy dependence of the 
s-wave potential. This leads only to a 10 % increase of Re Us at p = p^jl. 

B. Uncertainties from the Roper resonance 

There is another type of medium effect which has not been considered so far and will introduce additional un- 
certainties. This is related to the Roper excitation and its decay into nucleon and two pions in / = and s-wave. 
The Roper is the lightest resonance with the same quantum numbers as the nucleon and allows for a decay into a 
nucleon and two pions which are in isospin zero and s-wave relative to each other and also relative to the nucleon. In 
Ref. [83] the mechanism of Roper excitation from an isoscalar source and subsequent decay into two pions has been 
found dominant at low energies in the NNmr production for pions in / — 0. The isoscalar source can be 

described by an effective a exchange aNN* between the nucleons, whose strength has been fitted independently for 
the (a, a') reaction on a proton target [84]. Based on that finding, the relevance of this mechanism in nd scattering 
at low energies was also stressed in [40]. We shall also consider it here in connection with the s-wave pion- nucleus 
optical potential. 

For the present purposes the mechanism described above can be adapted by having the two pions, one in the initial 
state and the other one in the final state, as indicated in Fig. 26 on the left hand side. As the two pions are in a 
relative / = state, we obtain an isoscalar contribution to the ttN amplitude. The second nucleon line to which the 
isoscalar a couples is closed and gives a medium contribution to irN scattering. In the heavy baryon limit the diagram 
reduces to a point-like interaction of a pion with two nucleons as indicated in Fig. 26 on the right hand side. 

For the N*NTnr coupling an effective Lagrangian from Ref. [85] is used which leads to the effective vertex [83] 

-iSH^.j,^^ = -2i ^ U cl ^) (52) 

for 7r+7r+, tt^tt^ and tt^'tt'' and zero otherwise (note a minus sign in C2 with respect to [83] because now one of the 
pions is incoming). Here 101,102 are the energies of the pions. The values for the couplings are obtained in Ref. [83] 
from a fit to the experimental width of the N* decay into A^tt+tt" and Ntt^tt^ , c* = — 7.27GeV^^ and C2 = OGeV^^ 
For the N*aN coupling the effective vertex is —iAHa-NN" = iF{q)gaNN-' where g'^j^j^, /{Att) = 1.33 and F a form 
factor of the monopole type for the off-shell a with A^r = 1.7 GeV, rricr = 550 MeV. 

In the heavy baryon approximation we can put the external nucleons at rest and, thus, obtain for the elastic 
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FIG. 26: The Roper resonance in isoscalar, s-wave nN scattering in the medium. On the right hand side the interaction in the 
lieavy baryon limit is shown. 



scattering of a pion of any charge with a nucleon of any charge 



{-it) = 2 -2i -rf (cj - C2) {iFa{qa)gaNN') T i'i'P<y{lcr)9aNN) (Pp + Pn) (53) 

V J - TUN' -mi 

where the aNN couphng is the same as in the Bonn model [86] with g'^NN/i'^'^) " ^-^^ ^'^'^ Qo-nn and g^NN* have 
the same sign [84]. The contribution in Eq. (53) already contains the sum of the two diagrams on the left hand side 
of Fig. 26. The isoscalar modification is. thus, 

2 ■ 0.184 77V — . (54) 

2TTf^mi(mN - mN*)(m^ + rrifq) po 

At normal nuclear matter density p ~ pQ this leads to an isoscalar of Sbo = 0.188 to^^ which implies attraction. 
The result from Eq. (54) is huge compared to the isoscalar from the model of ttN interaction from Ref. [31] of 
be — — 0.0336m~^, see Tab. III. We can use the 6bo from Eq. (54) and calculate Us from Eq. (1). Then, already at 
tree level, one obtains the unrealistically large attraction of Re 115 ~ —2.5 ■ lO^MeV^. 

However, by turning the pion line around in the diagram of Fig. 26, we have implicitly changed the kinematics at 
which the above couplings, such as gcNN*, have been determined. The A^*(1440) is now off-shell by around 500 McV 
{E — ttin) which induces unknown theoretical errors in the calculation. 

Instead of the set (c| — — 7.27GeV~^, C2 = OGeV~^) one can consider the results from Ref. [33] which use the 
combination c| + C2 = (—1.56 ± 3.35) GeV~^ from Ref. [85] and then apply a resonance saturation hypothesis for 
the c* to be saturated by scalar meson exchange. Then the combination c| — C2 can be disentangled by the relation 
c*/c^ = 4.2/3.2 and our result with these values would change from ReHs = -2.5 • lO"' MeV^ to (-710 ± 1600) MeV^ 
at p = Po (compare to Figs. 24, 25). 

For p = po/2, the effective density felt by pionic atoms, the contribution is RcIIs = — 177± 400 MeV^, which, 
added to the results calculated before and shown in Fig. 25, leads to Rellg — 1030 ± 400 MeV^ for = + 10 
MeV (for the Coulomb shift) or ReHs = 1430 ± 400 MeV^ for k° ^ + 20 MeV. In both cases this band overlaps 
with the phenomenological fits, but the amount of theoretical uncertainty is indeed large. In order to account for all 
the uncertainties, we have taken the results for fc" = 16 MeV of Coulomb shift, corresponding to an average nucleus 
with Z = 56, A = 100, and have summed in quadrature the uncertainties from the Roper contribution and those of 
Figs. 17 and 19. This leads to the hatched band of theoretical values plotted in Fig. 27. We can see that the band 
is relatively broad, broader than for Imlls, since, as we have seen, there are more sources of uncertainty for ReHs. 
Altogether, with the realistic theoretical uncertainties accounted for, we find two bands for Re lis and Imlls which 
overlap with the empirical values needed to describe pionic atoms, for the effective density p ~ Pa/2 felt by pionic 
atoms. 



VIII. SUMMARY AND CONCLUSIONS 



The s-wave pion-nucleus optical potential has been calculated in a microscopical many-body approach by looking 
simultaneously at vertex and selfenergy corrections. We have, thus, taken a chiral unitarized rescattering approach 
that delivers a good description of vacuum data in the vicinity of the threshold and above. Subsequently, the medium 
corrections have been added to the vacuum model. Whereas Pauli blocking in rescattering generates repulsion, the 
pion polarization for intermediate pions, including ph, Ah and short-range correlations, is responsible for a moderate 
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FIG. 27: Uncertainties of the present study (hatched area). Phenomenological fits to pionic atom data as in Fig. 24. 



attraction. The model has been formulated for asymmetric nuclear matter, although vertex renormalizations and 
other corrections have only been evaluated for symmetric nuclear matter. 

For the Weinberg- Tomozawa term and the isoscalar contribution from the NLO chiral Lagrangian, in-medium vertex 
corrections, some of them novel, have been included. E.g., the Weinberg- Tomozawa term is increased, in agreement 
with recent analyses which include also data on deeply bound pionic atoms. 

We have also investigated vertex corrections for the NLO isoscalar nN amplitude and have found them relevant 
and of the same relative size as the renormalization of the isovector vertex. 

When these corrections are taken into account within the multiple scattering series of the Bethe Salpeter equation, 
together with selfenergy insertions in the intermediate states and Pauli blocking of the nucleons, we obtain an s-wave 
pion selfenergy, Xlg, in good agreement with empirical determinations for the the imaginary part and only qualitative 
for the real part. 

An important ingredient, already exploited in former works, has been the effect of the Coulomb shift in the argument 
of Xlg, which appears in the Klein Gordon equation. This effect leads to an increase of the repulsion in Re Us which 
brings the pion nucleus optical potential in better agreement with empirical determinations. 

We noticed that an increased bi in the medium led through rcscattering to an attraction in Re Ilg, in spite of the fact 
that the Pauli correction to these terms, given by the Ericson-Ericson formula, produces an increased repulsion. The 
rcscattering term with the non Pauli blocked part is larger than with the Pauli blocked one and of opposite sign, thus, 
altogether, rescattcring with an increased bi produces a net attraction. We also noted that the vacuum model for ttN 
has a certain freedom in the choice of parameters that induce changes in the real part of the rcscattering amplitude, 
but this is compensated by the NLO isoscalar terms in vacuum. We have tested that, with the simultaneous change 
in the medium of the isovector vertex and the NLO isoscalar term, the results for lis are stable with respect to these 
changes of the vacuum model. 

Another source of contribution to Re lis was a genuine term related to the Roper coupling to two isoscalar 
pions. This two nucleon term was determined, within uncertainties, in former studies of the -kcL interaction and we 
have used this information to evaluate the contribution to Re Hg, resulting in a band of values with relatively large 
uncertainties that one must accept. 

Altogether, we determine a band of results for Re Us which overlaps with empirical determinations in the region 
of interest of pionic atoms, p ^ po/2, and a much narrower band for Im Us which agrees with the also narrower band 
of empirical analyses in the same region. 

Along the work we have pointed out sources of uncertainties which we have quantified and summed in quadrature 
at the end. This leads to larger uncertainties in Re lis than were assumed in former studies, and that one must bear. 
Yet, within these admitted uncertainties, the results obtained represent a satisfactory description of the s-wave pion 
selfenergy for pionic atoms. 
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